Doppler Shift in Andreev Reflection from a Moving Superconducting Condensate in 

Nb/InAs Josephson Junctions 
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We study narrow ballistic Josephson weak links in a InAs quantum wells contacted by Nb elec- 
trodes and find a dramatic magnetic-field suppression of the Andreev reflection amplitude, which 
occurs even for in-plane field orientation with essentially no magnetic flux through the junction. 
Our observations demonstrate the presence of a Doppler shift in the energy of the Andreev levels, 
which results from diamagnetic screening currents in the hybrid Nb/InAs-banks. The data for con- 
ductance, excess and critical currents can be consistently explained in terms of the sample geometry 
and the McMillan energy, characterizing the transparency of the Nb/InAs- interface. 

PACS numbers: 74.45.-l-c, 74.50.-|-r, 74.81.-g, 85.25.Cp 



In recent years a detailed microscopic understand- 
ing of the proximity effect has emerged. There is now 
agreement that in highly transparent Josephson junc- 
tions formed by a metaUic weak hnk in good contact 
with two superconducting (SC) banks the supcrcurrent is 
carried by Andreev bound states (ABS)^. These come in 
pairs corresponding to the opposite directions of Cooper- 
pair transfer mediated by multiple Andreev reflection 
(MAR) at the SC / metal interfaces provided that the ac- 
quired quasi-particle phase is a multiple of 27r.— At cur- 
rents exceeding the critical current Ic{T,B) MAR be- 
tween the SC banks manifests itself in the current- voltage 
characteristics as subharmonic gap structures at volt- 
ages eVn = 2A/n, where A is the SC energy gap and 
n = 1,2, .... At higher voltage eV 3> 2A, I{V) becomes 
linear with an excess current /oxc = ^(^) ~ GnV deter- 
mined by a single Andreev reflection (AR) probability 
|a(e)p {Gn is the normal state conductance) 

Weak links formed by a two-dimensional electron gas 
(2DEG) in semiconductor quantum wells^ are of partic- 
ular interest because here the ballistic transport can be 
studied. In very high magnetic field perpendicular to the 
2DEG, theorj*^ and experiments^i^ have demonstrated 
Andreev transport via edge states. Indirect evidence for 
a strong magnetic-field effect on AR was experimentally 
found in antidot billiardsj^ The case of parallel field, with 
respect to the 2DEG, is equally intriguing: as ideally no 
magnetic flux threads the 2DEG, one may naively expect 
the Josephson current to survive up to the critical fields 
of the SC leads. This is not the case, but the under- 
lying mechanism of the supcrcurrent suppression is still 
unclear. This question, also relevant for other 2D hybrid 
systems,— is among the issues this study focuses on. 

In this Communication, Nb/InAs Josephson junctions 
of different width are studied in a four-terminal lead con- 
figuration within the 2DEG. This allows us to separately 
determine the transparencies of the InAs weak link and 
Nb/InAs interfaces, and identify an additional energy 
scale in the electronic spectrum of the hybrid SC termi- 
nals. We observe a very strong suppression of both the 



AR probability and supcrcurrent in weak magnetic fields 
of 4 and 100 mT for perpendicular and parallel orienta- 
tions, respectively. This unexpectedly rapid decay can be 
traced back to diamagnetic supercurrents in Nb, which 
break the time-reversal symmetry of Andreev reflection. 

The samples were prepared from a heterostructure con- 
taining an InAs-quantum well with thickness d^v = 15 nm 
and a mean free path around 3.7 fim, confined be- 
tween two AlGaSb layers. The electron density was 
7.8 • 10^^ m^^, resulting in a Fermi wavelength Xp = 
2TT/kp = 28nmii First, a constriction of width W with 
a four-terminal (4t) lead structure (shown in yellow in 
Fig. [T^) was patterned using electron beam lithograpy 
(EBL) and etched by reactive ion etching in a SiCl4- 
plasma. In a second EBL step two Nb stripes were de- 
posited onto the InAs (horizontal grey bars in Fig. [1^) 
after removal of the top AlGaSb layer and in-situ Ar-ion 
cleaning of the InAs surface. The Nb-stripes have a width 
Ws = 1 /ini and a thickness ^5 = 35 nm; their distance 
L = 600 nm defines the length L of the junction. This 
results in a Thouless energy of exh = hvp/L — 0.8 meV 
< A for all samples, vp being the Fermi velocity in the 
normal metal. We have prepared four samples on two 
chips in the same batch. The chips contained junctions 
with W = 0.5 and 1 /im (denoted as #1 and #2a), and 
W = 1 and 2/im (denoted #2b and #3), respectively. 

Figure [T]d shows the two-terminal (2t) resistance of 
samples 7^1, #2a, and #3 measured through the Nb 
stripes vs. temperature T. The resistance drops at the 
transition temperature of the Nb stripes around 8.3 K 
and again between 2 and 3K, where the constriction 
becomes superconducting. In Fig. [TJ: we plot the 2t- 
resistance measured across the Nb stripes and the 4t- 
resistance measured within the InAs layer at T = 10 K 
as a function of W. The 4t-resistance scales as 1/W, 
obeying the ballistic Landaucr-Biittikcr formula: 
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FIG. 1: (color online) (a) Scanning electron micrograph of 
sample #3 with width W = 2 fim and length L = 0.6 /xm. 
Regions with InAs two-dimensional electron gas are shown in 
yellow. Etched regions are light grey. The horizontal dark 
grey stripes are made from a 35 nm thick Nb. (b) R{T) for 
three samples with width W — 0.5, 1, 2 /im. (c) Resistance 
vs. 1/W at 10 K for the 4t-configuration with InAs leads (0) 
and the 2t-configuration with Nb leads (□). (d) Normalized 
supercurrent Ic vs. perpendicular magnetic field B± for dif- 
ferent widths W. The black curve is the standard Fraunhofer 
pattern fitting the data for sample #3 {W — 2 /im). 



slope of Gn{W) we extract (T) = 0.8. The 2t-resistance 
at 10 K is also proportional to 1/W, however, with an 
offset caused by the normal state resistance of the Nb- 
stripes and the Nb/InAs- interface resistance. Finally, we 
show in Fig. [T]i the interference pattern of the junctions 
in a transversal magnetic field. For the widest junction 
the critical current Ic{B) exhibits the expected Fraun- 
hofer patternfii For the narrower ones the higher order 
maxima are suppressed, and in the 0.5 /im-wide device 
Ic{B) resembles more a Gaussian rather than a Fraun- 
hofer pattern. This will be addressed below. 

Information about the spectral distribution of ABS is 
contained in the T-depcndence of the supercurrent. In 
Fig. [2] we present the critical current at i? = as a 
function of temperature for different widths W = 0.5, 1 
and 2 /xm of the InAs constriction. With decreasing tem- 
perature Ic{T) initially increases and then saturates at 
T < 1 K. The solid lines are fits based on the approach 
of Grajcar et alJ^ who adapted the scattering theory for 
the Josephson curreniji^ to 2D ballistic proximity struc- 
tures. Two parameters were used in these fits: the av- 
erage normal-state transmission probability (T) and the 
McMillan energy VmcM = tsn ■ hvp /'^dNr^ which is the 
Thoulcss energy of the InAs quantum well multiplied by 
the Nb/InAs interface transparency tsn- ^McM deter- 
mines the strength of the proximity-induced supercon- 
ducting correlations in the InAs regions underneath the 
Nb film, which act as effective superconducting termi- 
nals described by the Green's functions of the McMil- 




FIG. 2: Critical current vs. temperature T for different widths 
W together with best fits according to our model. Inset: Val- 
ues of (T) (0) and Tmcm (□) resulting from the fits. 



Ian model.— Clearly, the smaller tsn, the weaker the 
proximity effect. Because (T) influences mainly the sat- 
uration value of Ic(T — > 0), while TmcM determines the 
decay of Ic with T, these two parameters can be ex- 
tracted independently from the measured Ic{T). 

The fit parameters are given in the inset of Fig. [21 The 
values of (T) scatter by ±5 % around 0.9, which is slightly 
higher, but still close to (T) = 0.8 estimated indepen- 
dently from 4t-resistance ([T|) . The values of TmcM scatter 
by ±15 % around 0.9 meV. From this, we find t^s ~ 0.06 
which is much smaller than (T) . The high mobility of the 
InAs-quantum well provides the high transparency of the 
constriction, while the comparatively low Nb/InAs inter- 
face transparency is sufficient to transform the Nb/InAs- 
bilayer into an effective superconductor with the addi- 
tional energy scale TmcM (see also inset in Fig. [S]). 

Having all sample parameters fixed, we now turn to 
the main topic of our study, i.e. the magnetic field de- 
pendence of Andreev reflection. First, we examine IV- 
and dl/dF-characteristics. Figure [3] shows the differen- 
tial conductance dl/dV of sample #2a {W = 1 /.m) for 
different values of perpendicular field B±. At B± = 
subharmonic gap structures (see arrows) appear at inte- 
ger fractions Vn = 2A]\[i,/ne of 2A]s[b- The location of the 
steps vs. 1/n, plotted in the inset of Fig. [3l agrees very 
well with the value of Ajv6 = 1.35 mcV. The n = 1 step 
appears at lower voltage as predicted in Ref. [l6l . A tiny 
perpendicular field of 4mT washes out these features, 
vanishing completely at lOmT. Also, the enhancement 
of dl/dV at low bias is suppressed. It is known^iii that 
the subharmonic gap structure and low bias enhance- 
ment of dl/dV originate from multiple AR at finite V. 
The large number of steps independently confirms the 
high AR probability |ap expected from the high values 
of (T). The rapid smearing of these steps thus indicates 
a surprisingly strong magnetic field dependence of |ap. 

This observation is further substantiated by the behav- 
ior of the excess current lexc^ obtained by integrating 
dl/dV — Gn (see Fig. over voltage. Since /oxc results 
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FIG. 3: (color online) Measured differential conductance 
dl /dV of sample #2a {W = for different values of a per- 

pendicular magnetic field B± = 0, 4, 10, 40, and 100 mT, from 
top to bottom. Arrows indicate the subharmonic gap struc- 
ture induced by multiple Andreev reflections. Inset: Voltages 
of the subharmonic gap structures for all samples. The solid 
line indicates the expected slope for Ajvt, = 1.35 meV. 



from a single AR at V > 2Ani,, it is more robust than 
MAR whose amplitude is cx |ap". Figure [4] shows Icxc on 
a logarithmic scale for both field orientations. It is sup- 
pressed to ~ 50 % at Bj_ = 30 mT and B\\ = 300 mT. 
The field exceeds B± by an order of magnitude, but 
is still a factor of 5 smaller than needed for an appre- 
ciable reduction of the gap AjytiB) as follows from the 
pair-breaking theorjii^ (dashed line in Fig. 2]). 

We attribute the behavior of /oxc(-B) to the magnetic 
field suppression of the Andreev reflection in the InAs 
proximity regions underneatli the Nb films. In the pres- 
ence of the external field B = V x A(r), tire SC or- 
der parameter acquires an inhomogencous phase. 7(r) = 
— 27r$^^/'^ A(r')(ir', whose gradient ks = V7 induces 
diamagnetic currents in the entire Nb/InAs structure, 
including the InAs proximity layer (sec. Fig. |4|3 and c). 
If ki and k f denote the wave vectors of initial and final 
quasi-particle states in the AR-process, momentum con- 
servation requires k^ -I- kj = ks , since the Cooper-pairs 
absorbed by the moving SC condensate have momentum 
fi-kg. This leads to a Doppler shift ed = —fi^k±ks /2m,* 
of the energy of the AR quasi-particle (fcj^ is the trans- 
verse wave number) The inset in Fig. [5] illustrates the 
energy dependence of |a(e)p. The fact that each of the 
Nb/InAs-banks constitutes a single hybrid superconduct- 
ing terminal is illustrated the two BTK-like peaks in the 
spectum, one at A^vh and another at a smaller energy 
determined by Fmcm- Similar hybrid excitation spectra 
have recently been observed by scanning tunneling spec- 
troscopy in diffusive Al/Cu-structures^. The suppres- 
sion of |a(£)P occurs in the range |e| < ^McM, when ed 
becomes comparable to the McMillan energy FmcA/-— 

The solid lines in Fig. show /cxc (B) obtained from 
the BTK-model)^ in which we express the AR amplitude 
a(e) in terms of the McMillan's Green's function o^^i^^ 
accounting for the Doppler shift sd in the InAs proxim- 
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FIG. 4: (color online) (a) Measured excess current of sample 
#2a vs. perpendicular (o) and parallel (•) magnetic field. 
The solid lines result from our model including the Doppler 
shift eoiks) (see text) and the dashed line shows 7cxc due to 
the suppression of A]V6(-B||) only.— (b) and (c) One of the 
Nb/InAs terminals in parallel and perpendicular fields. Red 
arrows indicate the SC order-parameter phase gradient, kg = 
V7, with ks = —ir/W X <l?/"l?o determined by the magnetic 
fiux $ through the bilayer for given field orientation. 



ity layer. Using the parameters extracted from Fig.^ the 
model reproduces the measured curves quite well (except 
for the shoulder in /cxc(Sj_), which at present is not un- 
derstood). The model accounts for the ratio of the B^^ 
and B\\ scales. This has a simple interpretation if we 
notice that the relevant value of ks is related, via the 
circulation theorem, to the magnetic flux 4> threading 
the Nb/InAs bilayer for given fleld orientation (Fig. Hja 
and c): ks = —t:/W x $/$o- Neglecting the field in- 
homogeneity in Nb,— we find $ = B\\W{ds + (In) and 
$ = B±WWs for the parallel and perpendicular cases, 
respectively. Consequently, the expected field ratio is 
B\^/B± = Ws/{ds + d]^) ~ 20, which is within the mar- 
gins of the experimental uncertainty. 

The dramatic reduction of the AR amplitude a{B), 
inferred from the behavior of /cxc, can also explain the 
suppression of the higher lobes of the critical current in 
the Fraunhofcr pattern in weak perpendicular magnetic 
field (see Fig. [TJi), since the same a{B) determines the 
energies of the ABS and hence Ic{B)J^'^ We also note 
that the smaller the channel width W is, the higher are 
the values of B± required to reach $ = $0 in the junction 
and the stronger is the suppression oiIc{B) at this point. 

Finally, we examined the critical current Ic in a paral- 
lel magnetic field for samples #1 and #2a. Again, even 
for the perfectly parallel field there is the phase gradi- 
ent due to the vector potential (Fig. ^p) that induces 
a Doppler energy shift of the ABS. When it becomes 
comparable to FmcA/, counter-propagating ABS overlap 
in energy, and their contributions to Ic gradually cancel. 
Similar to /cxc, the suppression of Ic occurs at 100 mT. 
In Fig.[5]we plot Ic{B) for sample #1 for several almost 
parallel field orientations. The small perpendicular com- 
ponent of B still leads to a Fraunhofer-like minimum in 
/c, which allows a precise determination of the angle 
a between B and the junction plane. This can be de- 
scribed by h{B) = h{B\\) |sin(7r^Bi/$o)/(7r^5±/$o)|, 
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FIG. 5: Critical current vs. B for sample #1 at small angles 
a between B and 2DEG. The lines are theoretical fits (see 
text). The inset shows the double peak structure of |a(£)P 
characteristic of hybrid SC and its suppression by _B|| . 

where = B sin a, -B|| = B cos a and A is the junc- 
tion area. UnUke the standard Fraunhofer pattern^ii we 
take into account the dependence of the critical current 
Jc(i?||) on the parallel component of the field. This is 
done by straightforwardly generalizing the ballistic for- 
mula for the Josephson current in which we express the 
AR amplitude a in terms of the McMillan's Green's func- 
tions with the Doppler shift in the InAs proximity layers. 
The dependence Ic{B\\) makes the Fraunhofer-like oscil- 



lations decay much faster than the usual 1 /B± law, which 
is clearly seen in Fig. [5] For a < 1.3° no minimum is de- 
tected, while it is still observed for a wider sample #2a on 
the same chip (not shown). The curve for a « resulted 
from a careful maximization of Ic at fixed field, since 
in this case, evidently, no Fraunhofer minimum can be 
observed anymore. The fits (solid lines in Fig. [5|) are ob- 
tained by varying the effective thickness of the Nb/InAs 
bilayer, doff. We find des = + ds = "oO nm for sample 
#2a (not shown) in agreement with the geometry, while 
for sample #1 c?eff ~ 25 nm. The reason for this discrep- 
ancy is unclear, since the excess current data in Fig. 2] 
agree very well. The deviation of the data for a sa from 
the theoretical B±^ = curve is probably due to a weak 
perpendicular stray field of the Nb electrodes. 

In conclusion, we have studied the IV-characteristics 
and critical current of ballistic Nb/InAs Josephson junc- 
tions in parallel and perpendicular magnetic fields. The 
observed field dependence is much stronger than antici- 
pated from standard models and can be traced back to 
the Doppler shift of Andreev levels due to diamagnetic 
supercurrents in hybrid Nb/InAs contacts. Several dif- 
ferent aspects of the proximity effect are consistently and 
nearly quantitatively described by our theoretical model. 
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